On the uncertainty in the 0vfl(3 decay nuclear matrix elements 



m 
o 
o 

>>: 



(N 



> 

in 
o 
o 
in 
o 
m 
o 



o 

G 



X 



V. A. Rodin 1 , Amand Faessler 1 , F. Simkovic 2 , and Petr Vogel 1 ' 3,4 
1 Institute fur Theoretische Physik der Universitat Tubingen, D-72076 Tubingen, Germany 
2 Department of Nuclear Physics, Comenius University, Bratislava, Slovakia 
3 Department of Physics 106-38, California Institute of Technology, Pasadena, CA 91125, USA 
4 Institute of Particle and Nuclear Physics, Charles University, Prague, Czech Republic 

(February 8, 2008) 

The nuclear matrix elements M°" of the neutrinoless double beta decay (Of /3/3) are evaluated for 
76 Ge, 100 Mo, 130 Te, and 136 Xe within the Renormalized Quasiparticle Random Phase Approximation 
(RQRPA) and the simple QRPA. Three sets of single particle level schemes are used, ranging in size 
from 9 to 23 orbits. When the strength of the particle-particle interaction is adjusted so that the 
2f /3/3 decay rate is correctly reproduced, the resulting M°" values become essentially independent 
on the size of the basis, and on the form of different realistic nucleon-nucleon potentials. Thus, one 
of the main reasons for variability of the calculated M 0l/ within these methods is eliminated. 
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The discovery of oscillations of atmospheric [1] , solar 
[2] , and reactor [3] neutrinos shows that neutrinos have a 
non-vanishing rest-mass. However, since the study of os- 
cillations provides only information on the mass-squared 
differences, it cannot by itself determine the absolute val- 
ues of the masses (or even, at present, the sign of the 
mass-squared differences). Moreover, flavor oscillations 
are insensitive to the charge conjugation properties of 
the neutrinos, i.e. whether massive neutrinos are Dirac 
or Majorana particles. Study of the neutrinoless double 
beta decay (Of/3/3) is the best potential source of this cru- 
cial information on the Majorana nature of the neutrinos 
and on the absolute mass scale. 

The search for Of /3/3 decay has a long history. So far, 
the decay has not been seen, but impressive limits on its 
half-life have been established (for the review of the field, 
see e.g. [4] or [5]). With the progress in oscillation studies, 
which now established the mass scale m sca i e rs V Am 2 , a 
wave of enthusiasm emerged in the community to develop 
a new generation of experiments, sensitive to such mass 
scale. In that context it is crucial to develop in parallel 
methods capable of reliably evaluating the nuclear ma- 
trix elements, and realistically asses their uncertainties. 
The goal of the present work is to make a contribution 
to that endeavor. 

Clearly, the determination of the effective Majorana 
mass (m„) can be only as good as the knowledge of the 
nuclear matrix elements, since the half-life of the Of/3/3 
decay and (m„) are related by 
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where G 0v (Eq,Z) is the precisely calculable phase-space 
factor, and M 0l/ is the corresponding nuclear matrix el- 
ement. Thus, obviously, any uncertainty in M makes 
the value of (m u ) equally uncertain. In turn, 
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where on are the unknown Majorana phases. The ele- 
ments of the mixing matrix |£/ e i| 2 and the mass-squared 
differences Am 2 can be determined in the oscillation ex- 
periments. Using the present knowledge of these quanti- 
ties, (see, e.g., [6]) and limits on T]/ 2 , one could decide 
whether the neutrino mass pattern is degenerate, or fol- 
lows the inverse or normal hierarchies. If, on the other 
hand, the existence of the Of /3/3 decay is proven and the 
value of Ti/2 is found, a relatively narrow range of ab- 
solute neutrino mass scale can be determined, indepen- 
dently of the phases in most situations [4,6]. However, 
such an important insight is possible only if the nuclear 
matrix elements are accurately known. 

The nuclear matrix clement M 0ly is defined as 
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where (|/)) are the wave functions of the ground 
states of the initial (final) nuclei. The explicit forms of 
the operators M^,M^ T and M$ v are given in Ref. [7]. 
In comparison with most of previous Of/3/3 decay stud- 
ies [8-11] the higher order terms of the nucleon current 
are also included in the present calculation, resulting in a 
suppression of the nuclear matrix element by about 30 % 
[7]. Note that in the numerical calculation of M 0l/ here 
the closure approximation is avoided and the unquenched 
values gA = 1.25, gv = 1-0 are used. 

Two basic methods are used in the evaluation of M ", 
the quasiparticle random phase approximation (QRPA) 
with its various modifications and the nuclear shell model 
(NSM). These two approaches represent in some sense 
opposite extremes. 

In the QRPA one can include essentially unlimited 
set of single-particle states, but only a limited subset 
of configurations (iterations of the particle-hole, respec- 
tively two-quasiparticlc configurations). In the context of 
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QRPA several issues that have been raised, and deserve 
a systematic study: 

• For realistic interactions the QRPA solutions are 
near the point of the so-called collapse and thus its 
applicability is questionable. Numerous attempts 
have been made to extend the method's range of 
validity by partially avoiding the violation of the 
Pauli principle. Here we consider the simplest of 
them, the renormalized QRPA (RQRPA) [12,8]. 
We compare the results of RQRPA with those of 
the standard QRPA [9-11]. 

• The choice of the size of single-particle (s.p.) space 
is to some extent arbitrary, often dictated by con- 
venience. What effect this choice has on the M 0u 
values is the main thrust of the present work. 

• There are various forms of the nucleon-nucleon po- 
tential that lead to somewhat different forms of 
the resulting G-matrix . By comparing the re- 
sults obtained with three such potentials (Bonn- 
CD [16], Argonne [17], and Nijmegen [18]) we show 
that the resulting M 0v are essentially identical, and 
independent of the choice of the realistic nucleon- 
nucleon potential. 

In contrast, in the NSM one chooses a limited set of 
single-particle states in the vicinity of the Fermi level, 
and includes all (or most) configurations of the valence 
nucleons on these orbits in the evaluation of M 0l/ . The 
main open question in this approach is to determine the 
effects of the neglected single-particle states further away 
from the Fermi level. As shown below, we have also per- 
formed QRPA and RQRPA calculations with the set of 
s.p. states used usually in NSM. It appears that these 
methods, at least with the nucleon-nucleon interaction 
we used, are not applicable for such small s.p. bases. 

In what follows the Qvf3[3 decay nuclear matrix ele- 
ments M 0l/ for 76 Ge, 100 Mo, 130 Te, and 136 Xe are eval- 
uated. These nuclei are most often considered as candi- 
date sources for the next generation of the experimental 
search for Qvf3(3 decay. For each of them three choices 
of the s.p. basis are considered. The smallest one has 
9 levels (oscillator shells 7V=3,4) for 76 Ge, and 13 lev- 
els (oscillator shells iV=3,4 plus the / + h orbits from 
N = 5) for 100 Mo, 130 Te and 136 Xe. For the intermedi- 
ate size s.p. base the N = 2 shell in 76 Ge and 100 Mo are 
added, and for 130 Te and 136 Xe also the p orbits from 
TV = 5. Finally, the largest s.p. space [8] contains 21 lev- 
els for 76 Ge and 100 Mo (all states from shells N =1 — 5), 
and 23 levels for 130 Te and 136 Xe (N = 1 - 5 and i orbits 
from N = 6). Thus the smallest set corresponds to lhu> 



particle- hole excitations, and the largest to about 4hu> 
excitations. The s.p. energies have been calculated with 
the Coulomb corrected Woods-Saxon potential. 

It is well known that the residual interaction is an ef- 
fective interaction that depends on the size of the single- 
particle basis. Hence, when the basis is changed, the 
interaction should be modified as well. Here we propose 
a rather simple way to accomplish the needed renormal- 
ization. 

There are three important ingredients in QRPA and 
RQRPA. First, the pairing interaction has to be included 
by solving the corresponding gap equations. Within the 
BCS method the strength of the pairing interaction de- 
pends on the size of the s.p. basis. As usual, we multi- 
ply the pairing part of the interaction by a factor g pa i r 
whose magnitude is adjusted, for both protons and neu- 
trons, such that the pairing gap is correctly reproduced, 
separately for the initial and final nuclei. 

Second, QRPA equations of motion contain a block 
corresponding to the particle-hole interaction, renormal- 
ized by an overall strength parameter g p h- That param- 
eter is typically adjusted by requiring that the energy 
of some chosen collective state, often the giant GT reso- 
nance, is correctly reproduced. We find that the calcu- 
lated energy of the giant GT state is almost independent 
of the size of the s.p. basis and is well reproduced with 
g p h ~ 1. Accordingly, we use g p h = 1 throughout, with- 
out adjustment. 

Finally, QRPA equations of motion contain a block 
corresponding to the particle-particle interaction, renor- 
malized by an overall strength parameter g pp . (The im- 
portance of the particle-particle interaction for the f3 
strength was recognized first in Ref. [19], and for the 
(3(3 decay in [9].) It is well known that the decay rate 
for both modes of (3[3 decay depends sensitively on the 
value of g pp . In the following we use this property to find 
the value of g pp for each of the possible s.p. bases. The 
value of the parameter g pp is fixed in each case so that 
the known half-life of the 2vf3(3 decay is correctly repro- 
duced. The 2v(3f3 half-lives and average matrix elements 
collected in Table 1 of Ref. [4] are used, where the orig- 
inal references to the corresponding experiments can be 
found. A similar compilation of the 2vfij3 data can be 
found in Ref. [20]. The resulting adjusted values of g pp 
are shown in Table I for both the RQRPA and QRPA 
methods. One can see that as the basis increases, the 
effective g pp decreases 2 , as expected. 

The adjustment of g pp is a crucial point of the present 
work. Several studies of the sensitivity of the nuclear ma- 
trix elements M 0u to various modifications of the QRPA 
method as well as to the number of s.p. states and values 



1 Brueckner reaction matrix elements with these nucleon-nucleon potentials were calculated and used in Refs. [13-15]. 
2 For 136 Xe we use the lower limit on the 2v half-life for the adjustment of g pp ; hence further adjustment might be needed 
when the 2v half-life becomes known. 
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of g pp were made in the recent past [8,21,22]. Typically, 
these studies concluded that the values of M 0u vary sub- 
stantially depending on all of these things. Here we come 
to a different conclusion: by requiring that the known 
2v(3(3 decay half-life is correctly reproduced, and adjust- 
ing the parameter g pp accordingly, we remove much of the 
sensivity on the number of single-particle states, on the 
NN potential employed, and even on whether RQRPA 
or just simple QRPA methods are used. 

Clearly, the chosen procedure of finding the effective 
interaction is rather crude. Ideally, properly evaluated 
effective Hamiltonian, as well as the corresponding effec- 
tive (3(3 operator should be used in each case. However, 
as shown below, the chosen procedure appears to be suffi- 
cient to stabilize the values of the nuclear matrix elements 
M {)v . Note, that we did not use the quenched value of 
the axial current coupling constant gA, as is often done 
in studies of ordinary (3 decay. We are convinced that 



with the quenched gA the basic conclusion of our work 
will be similar, even though the values of g pp in Table I 
would be of course different. 

Having fixed the parameters of the effective Hamilto- 
nian we can proceed and evaluate the Qv(3(3 nuclear ma- 
trix elements M 0l/ and then the corresponding half-life 
(we list the phase-space factors in Table II in units of 
10~ 25 years for (ra„) in cV.) In evaluating M av we take 
the short-range correlations into account in the standard 
way, i.e. by multiplying the operators with the square of 
the correlation function [23] 

f(r) = 1 - e- ar \\ - br 2 ); a = 1.1 fnT 2 , b = 0.68 fnT 2 . 

(4) 

Note that the effect of short-range correlations reduces 
the matrix elements M 0u by a factor of about two, in 
agreement with other evaluations. 



TABLE I. Values of the effective particle-particle strength parameter g pp for the different nuclei and different sizes of the 
single-particle space. The values in col. 3 correspond to the minimal s.p. space, in column 4 to the intermediate one, and in 
col. 5 to the largest s.p. space. In each of these there are three entries corresponding to the G-matrix based on the Bonn, 
Argonne, and Nijmegen potentials (in that order). For every considered nucleus the upper line was obtained in RQRPA and 
the lower one in the simple QRPA. In column 2 we give the corresponding half-life of the 2u(3(3 decay used in the adjustment 
(see text, for 136 Xe the lower limit is used). 



Nucleus 


7\ 2 ; 2 (in 10 20 y) 


min. s.p. space 


intcrm. s.p. space 


largest s.p. space 


v<i Gc 


13.0 


0.99 1.12 1.07 


0.88 1.00 0.95 


0.79 0.88 0.84 






0.85 0.96 0.91 


0.75 0.85 0.81 


0.65 0.72 0.69 


100 Mo 


0.08 


1.21 1.35 1.30 


1.09 1.21 1.17 


1.00 1.10 1.07 






0.94 1.03 1.00 


0.85 0.93 0.90 


0.78 0.85 0.83 


i a o Te 


27.0 


0.97 1.10 1.05 


0.90 1.01 0.97 


0.84 0.94 0.90 






0.84 0.94 0.90 


0.78 0.86 0.83 


0.71 0.78 0.76 


136 Xe 


8.1 


0.82 0.93 0.89 


0.77 0.87 0.83 


0.72 0.82 0.78 






0.73 0.83 0.79 


0.68 0.77 0.74 


0.62 0.70 0.67 



TABLE II. Values of the nuclear matrix element |M°"| for the different nuclei and different sizes of the single-particle space. 
In column 2 is corresponding phase space factor G 0l/ . For the explanation of the notation in columns 3-5 see the caption of 
Table I. 



Nucleus 


G°" (in 10 _25 y-eV~ 2 


min. s.p. space 


intcrm. s.p. space 


largest s.p. space 


76 Go 


0.30 


2.41 2.37 2.35 
2.68 2.62 2.65 


2.52 2.44 2.47 
2.81 2.74 2.72 


2.32 2.34 2.35 
2.62 2.64 2.64 


100 Mo 


2.19 


1.08 1.08 1.05 
1.19 1.22 1.19 


1.12 1.14 1.08 
1.25 1.28 1.20 


1.28 1.34 1.27 
1.38 1.41 1.40 


iao Te 


2.12 


1.42 1.32 1.34 
1.50 1.42 1.40 


1.40 1.33 1.32 
1.46 1.42 1.39 


1.17 1.13 1.14 
1.22 1.17 1.15 


i3B Xe 


2.26 


1.05 0.99 0.99 
1.11 1.02 1.05 


1.09 1.03 1.04 
1.15 1.06 1.06 


0.92 0.83 0.86 
0.98 0.90 0.91 



As pointed out earlier, in the nuclear shell model an 
even smaller set of single-particle states is used corre- 
sponding to Ohui. This choice reflects the practical com- 
putational limitations in handling the extremely large 



number of possible configurations, while it seems to be 
sufficient to describe the spectroscopy of low-lying nu- 
clear states. In the NSM evaluation of the (3(3 decay 
rates [24] four s.p. orbits (/ 5 /2,P3/2,f>i/2, 59/2) were used 
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for 76 Ge and a five orbits (c? 5 / 2 , ^3/27 s i/2j 57/2; ^11/2) f° r 
130 Te and 136 Xe. These s.p. sets are free of the spurious 
center-of-mass states, but obviously miss a large part of 
the GT strength as well as of the strength corresponding 
to the higher multipoles. In order to describe GT transi- 
tions between low-lying states in the NSM, it is necessary 
to quench the corresponding strength. This is most con- 
veniently formally achieved by using g& = 1.0 instead of 
the free nucleon value of gA = 1-25. We follow this pre- 
scription in our attempt to use this smallest s.p. space, 
and only there. 

It appears that it is impossible to describe the 2v(3(3 
decay in such s.p. space using QRPA or RQRPA, and the 
nucleon-nucleon potentials employed in this work. One 
would have to renormalize the particle-particle block too 
much, with g pp ~ 2.0, unlike the rather modest renor- 
malization shown in Table I. With such large value of 
g pp the interaction is too far removed from the G-matrix 
used in the rest of this work. Therefore, one cannot ex- 
pect to obtain sensible (V matrix elements. In fact, we 
obtained very small matrix elements in this case for 130 Te 
and 136 Xe, while, perhaps accidentally, for 76 Ge they are 
in a crude agreement with the NSM result [24]. 
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FIG. 1. Dependence of the matrix elements MivQQ (left 
scale, dashed lines) and Mo v pp (right scale, full lines) on the 
parameter g pp . Calculations were performed for 9 and 21 
s.p. levels for 76 Ge as indicated; the Nijmegen potential and 
RQRPA method were used. The thin dotted horizontal line 
indicates that by fixing g vp to reproduce the experimental 
value M 2 „f3f3 = 0.15 MeV - the value of M 0v gg is also stabi- 
lized. 



We list the results with the three larger single-particle 
bases in Table II which represents the most significant 
part of the present work. As one can see by inspecting 
the entries, one can draw two important conclusions: 

• The resulting M " do not depend noticeably on the 
form of the nucleon-nucleon potential used. That 
is not an unexpected result. 

• Even more importantly, with our choice of g pp the 
results are also essentially independent on the size 
of the s.p. basis. This is a much less obvious and 
rather pleasing conclusion. It can be contrasted 
with the result one would get for a constant g pp 
independent on the size of the s.p. basis. The val- 
ues of M 0u differ then between the small and large 
bases by a factor of two or more. 

The effect of the g pp adjustment is illustrated in Fig. 
1, showing that our procedure leads to almost constant 
M 0v matrix elements. On the other hand, by choosing a 
fixed value of g pp the resulting M " matrix elements for 
9 and 21 s.p. levels would differ substantially. 

The entries in Table II are relatively close to each 
other. To emphasize this feature, each calculated value is 
treated as an independent determination and for each nu- 
cleus the corresponding average (M°") matrix elements 
(averaged over the three potentials and the three choices 
of the s.p. space) is evaluated, as well as its variance a 
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{M Uv )Y, (N = 9). (5) 



These quantities (with the value of a in paretheses) are 
shown in Table III. Not only is the variance substantially 
less than the average value, but the results of QRPA, al- 
beit slightly larger, are quite close to the RQRPA values. 
The averaged nuclear matrix elements for both methods 
and their variance are shown in Fig. 2. 

Combining the average {M "} with the phase-space 
factors listed in Table II the expected half-lives (for 
RQRPA and (m„) = 50 mcV, the scale of neutrino masses 
suggested by oscillation experiments) are also shown in 
Table III. These predicted half-lives are a bit longer (par- 
ticularly for the last three nuclei on our list) then vari- 
ous QRPA calculations usually predict. They are faster, 
however, then the shell model results of Ref. [24]. 

TABLE III. Averaged Ov/3/3 nuclear matrix elements 
(M l,I/ ) and their variance a (in parentheses) evaluated in the 
RQRPA and QRPA. In column 4 the 0u(3f3 half-lives evalu- 
ated with the RQRPA average nuclear matrix element and 
for the (m v ) = 50 meV are shown. 



Nucleus 


RQRPA 


QRPA 


T 1/2 (in 10 27 y for (m„) = 50 mcV 


76 Go 


2.40(0.07) 


2.68(0.06) 


2.3 


100 Mo 


1.16(0.11) 


1.28(0.09) 


1.4 


l3U Tc 


1.29(0.11) 


1.35(0.13) 


1.1 


136 Xc 


0.98(0.09) 


1.03(0.08) 


1.9 
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• RQRPA 
O QRPA 



Ge Mo Te Xe 

FIG. 2. Average nuclear matrix elements (M°") and their 
variance for both methods and for the four considered nuclei. 



Given the average nuclear matrix elements in Table 
III and the phase space factors in Table II one can find 
a limit (or actual value) of the effective neutrino mass 
(m„) from any limit (or value) of T1/2 from 



(cV) = (T 1/2 x G 



-1/2 
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(6) 



In conclusion, we have developed a "practical" way of 
stabilizing the values of the Qvf3f3 nuclear matrix elements 
against their variation caused by the modification of the 
nucleon-nucleon interaction potential and the chosen size 
of the single-particle space. We have also shown that 
the procedure yields very similar matrix elements for the 
QRPA and RQRPA variants of the basic method. Even 
though we cannot guarrantee that this basic method is 
trustworthy, we have eliminated, or at least greatly re- 
duced, the arbitrariness commonly present in the pub- 
lished calculations. 
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